Geometric Time Delay Interferometry 
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The space-based gravitational-wave observatory LISA, a NASA-ESA mission to be launched after 
2012, will achieve its optimal sensitivity using Time Delay Interferometry (TDI), a LISA-specific 
technique needed to cancel the otherwise overwhelming laser noise in the inter-spacecraft phase 
measurements. The TDI observables of the Michelson and Sagnac types have been interpreted 
physically as the virtual measurements of a synthesized interferometer. In this paper, I present 
Geometric TDI, a new and intuitive approach to extend this interpretation to all TDI observables. 
Unlike the standard algebraic formalism, Geometric TDI provides a combinatorial algorithm to ex- 
plore exhaustively the space of second- generation TDI observables (i.e., those that cancel laser noise 
in LIS A- like interferometers with time-dependent armlengths). Using this algorithm, I survey the 
space of second-generation TDI observables of length (i.e., number of component phase measure- 
ments) up to 24, and I identify alternative, improved forms of the standard second-generation TDI 
observables. The alternative forms have improved high-frequency gravitational- wave sensitivity in 
realistic noise conditions (because they have fewer nulls in the gravitational-wave and noise response 
functions), and are less susceptible to instrumental gaps and glitches (because their component phase 
measurements span shorter time periods). 

PACS numbers: 04.80.Nn, 07.60.Ly, 95.55.Ym 



I. INTRODUCTION 

The Laser Interferometer Space Antenna (LISA) is a 
joint NASA-ESA deep-space mission to be launched af- 
ter 2012, aimed at detecting and studying gravitational 
waves (GWs) with frequencies between 10" 5 and 1CU 1 Hz 
. LISA will provide access to GW sources that are out- 
side the reach of ground-based intcrfcromctric GW detec- 
tors 0, such as the binaries of compact stellar objects 
in our galaxy, the mergers of massive and supermassive 
black holes, and the gravitational captures of compact 
objects by the supermassive black holes at the center of 
galaxies 

LISA consists of three widely separated spacecraft, 
flying around the Sun in a quasi-equilateral triangular 
configuration and exchanging phase-coherent laser sig- 
nals. LISA relies on picometer interferometry to mea- 
sure GWs as modulations in the distance between the 
spacecraft. The greatest challenge to achieving this mea- 
surement is the phase noise of the LISA lasers, which 
is larger than the GW-induced response by many orders 
of magnitude, and which cannot be removed by conven- 
tional phase-matching interferometry because the LISA 
armlengths are grossly unequal, and changing continu- 
ously. Time Delay Interferometry (TDI), developed by 
J. W. ArmstrongjF. B. Estabrook, M. Tinto, and others 
H H H HH |1 [IS El, is the LISA-specific technique 
that will be used to combine the laser-noise-laden one- 
way phase measurements performed between the three 
spacecraft 1 into virtual interferometric observables where 



laser noise is reduced by several orders of magnitude. 

TDI was initially developed using ad hoc algebraic rea- 
soning for the case of a stationary LISA configuration 
with unequal but constant armlengths {first- generation 
TDI, see 0,0). It was later modified to work also in 
the case of a rotating LISA constellation (modified TDI, 
see 0, U 0, EI) and of linearly changing armlengths 
(second- generation TDI, see |lCt 111)). First-generation 
and modified TDI were given a rigorous mathematical 
foundation in the theory of algebraic syzygies on mod- 
uli 0, providing tools to generate all possible TDI ob- 
servables, and to determine which observables are op- 
timally sensitive to GWs 0. Unfortunately, this al- 
gebraic treatment cannot be extended easily to second- 
generation TDI, which is the version that must be used 
in practice. 

In this paper, I give a new derivation of first- 
generation, modified, and second-generation TDI, using 
a geometric approach that emphasizes the physical inter- 
pretation of TDI observables as synthesized interferomet- 
ric measurements 0, E3, El , extending it to all known 
observables. What is more, this geometric approach to 
TDI (in short, Geometric TDI) allows the exhaustive 
enumeration of all TDI observables of any length, and 
it leads to alternative, improved forms of the standard 
TDI observables, characterized by better GW sensitivity 
at high frequencies in realistic noise conditions, by lesser 
demands on the measurement system, and by reduced 
susceptibility to gaps and glitches. 

More specifically, all TDI observables display nulls in 



A variant of the technique uses combinations of one-way and two- 
way phase measurements, generated by locking five of the six 
LISA lasers to the last one, as described by Tinto and colleagues 



|12| . In this paper we shall consider only the one-way formalism, 
but our results could be applied with superficial modifications 
also to the two-way variant. 
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FIG. 1: In this idealization of the basic time-transport observ- 
able used with LISA, two ideal clocks travel along geodesies, 
with clock 1 continuously transferring its proper time to clock 
2 by way of pulsed light signals. GWs are measured as the 
fluctuations in the time of flight between the clocks (see main 
text). 



Geometric TDI: in Sec. Ill Al I introduce the basic GW- 
sensitive phase measurement; in Sec. Ill Bl I discuss its 
integration into laser-noise-canceling observables accord- 
ing to the Geometric TDI principle; in Sees. HTCl and 
III Dl I give a new derivation of the observables of first- 
generation, modified, and second-generation TDI, and I 
interpret them geometrically; in Sec. Ill El I show how to 
enumerate exhaustively all possible observables by repre- 
senting them as link strings; last, in Sec. Ill F| l extend our 
formalism, developed for simplicity by considering only 
three independent LISA lasers, to the realistic case of six 
lasers. Section ILTT1 reports on the exhaustive survey of all 
second-generation TDI observables consisting of up to 24 
separate phase measurements: in Sees. MI Al and MI Bl 
I discuss the alternative forms of the standard second- 
generation TDI observables, and present their practical 
advantages for the implementation of TDI; in Sec. MI d I 
describe the previously unknown second-generation TDI 
observables of length 18 and more. Last, Sec. II VI presents 
my conclusions. The appendices contain rules and proofs 
omitted from the main text, and explicit algebraic expres- 
sion for the second-generation TDI observables of length 
16. 

As customary, I set G — c — 1 except where specified 
otherwise. 



their noise and GW responses at frequency multiples of 
the inverse arm-crossing light times. 2 Because these zeros 
occur at the same frequencies and with the same orders 
for noise and GWs, the ideal GW sensitivity after suc- 
cessful laser-noise suppression is finite, and comparable 
to the sensitivity at nearby frequencies. The actual sen- 
sitivity, however, is likely to be degraded , ei ther because 
noise leaks into the nulls from the sides |15|. or because 
the measurement system has insufficient dynamical range 
to resolve the tiny signals within the nulls. This prob- 
lem is mitigated with the alternative observables, which 
have half as many response- function nulls as the standard 
forms. 

In addition, because the alternative observables are, 
as it were, folded versions of their standard forms, they 
have a smaller temporal footprint: that is, they are writ- 
ten as sums of one-way phase measurements that span 
a shorter time period. This property can be advanta- 
geous in the presence of instrumental gaps or glitches, 
which would then contaminate a smaller portion of the 
data set; a reduced temporal footprint means also that a 
shorter continuous set of phase data needs to be collected 
before TDI observables can begin to be assembled. 

This paper is organized as follows. SectionMldescribes 



2 The responses are exactly null only in the limit of equal LISA 
armlengths. For realistic, time-evolving LISA geometries the 
nulls are spread into narrow dips; however, these are deep enough 
that the qualitative discussion to follow still applies. 



II. A GEOMETRIC VIEW OF TIME-DELAY 
INTERFEROMETRY 

How is LISA an interferometer other than by name? 
The loosest dictionary definition of "interferometer" 
(something like "a device that combines the signals ra- 
diating from a common source, and received at different 
locations, or at the same location after traveling differ- 
ent paths" ) does not seem to apply to LISA, whose TDI 
GW observables are combinations of the phase-difference 
measurements between as many as six laser sources. In 
fact, intcrfcromctry is not needed, strictly speaking, to 
measure GWs, but only to remove the otherwise deafen- 
ing phase noise produced by the LISA lasers. The basic 
principle of GW measurement employed by LISA is non- 
interferometric, as we can see from the idealized experi- 
mental setup of a time-transport link between two ideal 
clocks (see Fig. P). 



A. The basic time-transport observable 

Consider a plane GW propagating across the 
Minkowski background geometry, and written in the 
transverse-traceless gauge 0. The wave is traveling 
along the x direction, and has "+" polarization along 
the y and z directions. We can then write the spacetime 
metric as r)p V + /iJJ, where 

h^y =h + (t + x)[e zz -e yy ]. (1) 
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Consider also two ideal clocks 1 and 2, marking their 
proper times t\ and £ 2 ? an d sitting at constant spatial 
coordinates pi = {0,0,0} and p 2 = {0,0, L} in the TT 
frame. In this gauge, constant-coordinate worldlines are 
geodesies, so the effect of the GWs is not to exert forces 
(as it were) on test particles, but to modulate the dis- 
tance between them. By way of light signals, clock 1 is 
continuously sending its time t\ to clock 2, where t\ is 
compared with the local time t 2 , yielding the difference 

At 12 = tj,(t) - ti(t - L 12 (t)) = L 12 (t). (2) 

Here t is the TT coordinate time and Li 2 (t) is the time of 
flight between the two clocks, as experienced by the laser 
pulse that arrives at clock 2 at time t. We are assuming 
that the two clocks have been synchronized so that in the 
absence of GWs they both mark the coordinate time t. 
To first order in h +1 L\ 2 {t) is 



L 12 (t)=L + 



h + (t)dt; 



(3) 



t-L 



the x coordinate dependence of the GW does not ap- 
pear in Eq. © because the two clocks sit on the same 
constant-x wavefronts. If the rates of the two clocks 
remain synchronized (dti/dt — dt 2 /dt), then the time 
derivative dAti 2 /dt is directly proportional to difference 
of the GW strains at the events of pulse reception and 
emission, 

This is our GW observable. In the Fourier domain, 
the (power) response function of dAti 2 /dt to GWs is 
exp(-27ri/L)| 2 /4 = sin 2 (tt/£) . Thus dAt 12 /dt is 
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insensitive to GWs of frequencies / <SC 1/L or / ~ k/L 
(with integer k). Expression 0} is the basic building 
block used to derive the LISA response to GW waves, as 
well as the Doppler response used in spacecraft-tracking 
GW searches [l7| , and the timing-residual response used 
in pulsar-timing searches |l8j . 

To relate this idealized experimental setup to LISA, 
we replace the ideal clocks with the LISA lasers, and ob- 
tain proper time by dividing the lasers' phase by their 
frequency. Each of the three LISA spacecraft contains 
two optical benches oriented facing the other two space- 
craft; on each bench, the appropriately named phaseme- 
ters compare the phase of the incoming lasers against 
the local reference laser. As written, Eq. (@J involves a 
comparison of laser frequencies: we choose to develop our 
arguments in terms of these, since it is more convenient 
to deal with instrumental responses that are directly pro- 
portional to the physical observable of interest (the GW 
strain) rather than to its time integral. Generalizing 
Eq. 0} to arbitrary plane-GW and spacecraft geometries, 
and adopting a LISA-specific language, we come to the 
Estabrook-Wahlquist two-pulse response 0] 

2 l-n^(i)fc 

m 

(5) 



with this indexing, Eq. (0 describes the frequency- 
difference measurement performed on spacecraft 2 to 
compare the local laser to the laser incoming from space- 
craft 1. In this equation: 

• kj is the spatial propagation vector of the plane 
GW; 

• p\ (t) and p 2 (t) are the spatial TT coordinates of 
the two spacecraft; 

• t is the time of pulse reception, and therefore of 
measurement; 

• i s is the time of pulse emission, as determined im- 
plicitly by \p s 2 (t) - p\{t s )\ =t-t B ; 

• i s the unit vector along the trajectory of 
the light pulse (labeled by the time of reception t) , 
given by nf 2 (i) = (p s 2 (t) - pf(i.))/|p5(t) - pf(i.)|. 

Equation 10 is known as the two-pulse response because 
an impulsive GW is registered twice in each (t) observ- 
able, once when it impinges on the emitting spacecraft i, 
and once, a time L later, when it impinges on the receiv- 
ing spacecraft j. 

In the literature on TDI, it is customary to label the 
LISA arms by the index of the opposite spacecraft. We 
shall do so in this paper, using primed and unprimed in- 
dices to denote the oriented LISA arms (with orientation 
following the direction of laser transmission) according 
to the convention 3 {1,2,3} = {3 -> 2, 1 -> 3, 2 -> 1} 
and {l',2',3'} = {2 3, 3 1, 1 -> 2}. We shall 
then denote by L/ the propagation time experienced by 
a laser pulse traveling along the oriented arm I. We shall 
also find it useful, at times, to augment the phase- 
measurement notation with a middle index, correspond- 
ing to the oriented arm traversed by the laser pulse being 
measured. (In fact, the primed or unprimed middle in- 
dex would be sufficient to identify the phase measurement 
completely, and we shall exploit this property in Sec. Ill El 
when we represent TDI observables as link strings.) See 
Fig. [5] for an example of this convention at work. 



B. The Geometric TDI principle 

Unfortunately, GWs cannot be read off directly from 
the yij measurements, because the fluctuations Cj(t) of 
the laser frequencies (i.e., the laser phase noises) come 
into the Uij(t) as 



y ij (t) = C i (t B )-C j (t) + GWB; 



(6) 



3 See the TDI Rosetta Stone, first printed in Ref. Il'il and updated 
at www . vallis . org/tdi for a mapping between the different con- 
ventions used in the literature. 
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FIG. 2: These time-delayed sums and differences of two yij 
measurements cancel laser phase noise at time t. In all of 
them, two laser pulses arrive at, or depart from, the same 
spacecraft at time t. 



the LISA lasers have Ci(t) of (single-sided, square-root) 
spectral density ~ 30Hz/\/Hz, several orders of magni- 
tude stronger than the weakest GWs detectable by LISA, 
which are at the level of the other two fundamental LISA 
noises (known together as secondary noises): the shot 
noise at the phasemeter, as determined by the power of 
the lasers and by the distance between the spacecraft, 
and the acceleration noise of the proof masses enclosed 
within each optical bench, which are used to reference 
the frequency measurements to freely falling worldlines. 
[Equation JBJ assumes that a single laser is being used on 
each spacecraft; it is pedagogical to consider this simpli- 
fied case first, but we shall generalize our discussion to 
the realistic case of six LISA lasers in Sec. Ill Fl ] 

Canceling laser phase noise is where interferometry 
comes to the rescue. Look at Fig. |2 for combinations of 
yij measurements in which two laser pulses arrive simul- 
taneously at spacecraft 1 at time t, depart simultaneously 
from spacecraft 1 at time t, or arrive and depart simulta- 
neously to and from spacecraft 1 at time t. We subtract 
the y^ measurements, represented graphically by arrows, 
when they share the same event of emission or reception 
(i.e., when their arrowtails or arrowheads meet), and we 
add them when the receiving spacecraft of one measure- 
ment is the emitting spacecraft of the other (i.e., when 
arrowtail follows arrowhead). In all of these combina- 
tions, the laser-frequency noise C\(t) generated at time 
t on spacecraft 1 is canceled out by entering twice with 
opposite signs; however, GWs are not canceled (not even 
at time t), because they come into Eq. © with differ- 
ent n™-dependent projection factors. The combinations 
of Fig. [2] do still contain frequency noise from lasers 2 
and 3, and from times other than t; it is however a sim- 
ple leap to cancel even those by arranging together more 
measurements. We then formulate a 

Geometric TDI principle: to obtain a laser-noise- 



canceling GW observable, line up arrows (i.e., 2/y mea- 
surements) head to head, tail to tail, or head to tail, 
creating a closed loop that cancels laser noise at all pulse 
emission and reception events. If no arrowhead or ar- 
rowtail is left unpaired, the closed loop represents a lin- 
ear combination of delayed y^ measurements that com- 
pletely cancel the three laser noises C,*(t). 

Remarkably, it is usually possible to interpret each 
closed-loop combination as the interferometric measure- 
ment performed by comparing the phases of laser beams 
that follow the paths marked by the arrows. Let us see 
an example. The arrows of Fig. (left panel) reproduce 
the paths followed by light in an equal-arm Michelson 
interferometer; operating in analogy with Fig. [21 and at- 
tributing the time t to the final common event of recep- 
tion at spacecraft 1, we write the corresponding algebraic 
expression 

yi 3 > 2 {t ~ L 3 ) + J/231 (*) - 2/32'i(*) ~ Vi23{t - L 2 >). (7) 

Here the sum of two time-consecutive 2/ij observables, 
such as y\ 3 ' 2 {t — L3) and 2/231 (i) [here L 3 is the light- 
travel time between spacecraft 2 and 1], simulates the re- 
flection of the laser off a mirror: in terms of laser phases, 
we see that the integral of this sum reproduces the to- 
tal phase shift accumulated along the path I — > 2 — > 1. 
By contrast, the head-to-head difference of two such dou- 
ble arrows simulates a photodetector: it reproduces the 
difference of the phase shifts accumulated along the two 
paths. All in all, Eq. J7J shows that the combination 
of four (one-way) 2/y measurements can synthesize the 
phase-difference output of a Michelson interferometer, 
as emphasized by Tinto and Armstrong Q, and shown 
graphically by Shaddock ^(j an d Summers 0[. Insert- 
ing the laser noises C,(i) in Eq. J7J), we get 



[C^t-Ls-Lv)- C 2 (t-L 3 )] 

+ [C 2 (t-L 3 ) -Ci(t)] 
- [C 3 (t -!*,)-&$)] 
\C x (t - L 2 , - L 2 ) - C 3 {t-L 2 ,% 



(8) 



which sums up to zero in interferometer geometries 
where 4 L 3 + L 3 > = L 2 i + L 2 : our equal- arm Michelson 
combination is then truly laser-noise canceling. It is how- 
ever sensitive to GWs, as can be seen by inserting Eq. 
© in Eq. J7J). 

More generally, we can set three simple rules to turn a 
closed arrow loop into a combination of y^ measurements 
that cancels laser noise: 

1. start at any spacecraft, and write down the appro- 
priate y^ for each arrow, following the loop (going 



4 If the interferometer is rotating, the Sagnac effect |21| introduces 
a distinction between the light-travel times L; and Ly in the two 



directions 8, 9]. 
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Equal-arm Michelson Unequal-arm Michelson 




FIG. 3: Left. — The arrows of this closed loop reproduce the 
paths of light in a standard equal-arm Michelson interferom- 
eter, and the corresponding time-ordered sum of phase mea- 
surements [Eq. 0] reproduces the phase-difference output of 
the interferometer. Right. — For unequal armlengths, laser- 
phase-noise cancellation can be recovered by having both in- 
terfering beams travel along each arm once, building up the 
same light-travel time. Compare with Fig. of Ref. |ic||. 



along or against the direction of each arrow) until 
all arrows are used up (if more than two heads or 
tails meet at any spacecraft, different visiting or- 
ders will yield different observables); 

2. use a plus (minus) sign for arrows followed along 
(against) their direction; 

3. give time arguments to the y^ , remembering that 
measurements are always made at the receiving 
spacecraft (at the arrowhead) ; use the nominal time 
t for the first y^j, and then add (subtract) the ap- 
propriate Li for each arrow followed along (against) 
its direction. 



C. The observables of first-generation TDI 

The first laser-noise-canceling combinations for LISA 
were discovered using an algebraic (rather than geomet- 
ric) approach, matching up delayed yij measurements 
in such a way that all laser-noise terms would cancel. 
Usingthis procedure, Tinto, Armstrong, and Estabrook 
[USB obtained expressions for first-generation TDI ob- 
servables, which cancel laser noise in static unequal-arm 
geometries. These observables are sums of either six or 
eight delayed j/y measurements (for short, links). See 
Fig. m 

The 6-link observables a, (3, 7 (mapped into each other 
by relabeling the spacecraft cyclically) use all six LISA 
oriented arms, and measure the phase difference accu- 
mulated by two laser beams traveling around the LISA 
array in clockwise and counterclockwise directions: thus, 
they behave much like a Sagnac interferometer, and are 
known as Sagnac observables. A related 6-link combi- 
nation, the symmetrized Sagnac observable £, has the 
useful property of being relatively insensitive to GWs in 



the low- frequency limit. 5 

The 8-link observables X, Y, Z (also mapped into 
each other by cyclic spacecraft relabelings) use two of 
the LISA arms in the two directions. They are unequal- 
arm generalizations of the Michelson observable of Eq. 
|Q: for unequal arms, the latter would fail to cancel the 
laser-noise terms from the tails of the two paths, because 
L3+L3' ^ L-2'+L2- The solution is to have both paths go 
through each arm once (hence the eight terms), building 
up the same light-travel time (see the right panel of Fig. 
0. Related 8-link combinations, known as observables of 
the U, P, and E type, use different sets of four oriented 
arms out of six, and have GW sensitivity comparable to 
the Michelson combinations [5|,|6j. 

Prior to my work, it was unclear whether the P-type 
and P-type observables could be interpreted as synthe- 
sized interferometric observables. 6 In Fig. 0] we show 
that this is possible if we identify jour distinct laser 
beams, paired in alternative ways to cancel laser noise 
at the path tails (dots) and path heads (ending arrows). 
The two path origins are not simultaneous, and neither 
are the two path endings. 

The symmetrized Sagnac observable £, which also de- 
fies explanation as a two-beam synthesized interferom- 
eter, can be interpreted clS el SIX— beam interferometer, 
whereby two different pairings explain the cancellation 
of laser noise at emission (dots) and reception (arrows). 
Yet another pairing, shown by the thin diagonal lines in 
Fig. explains why £ is relatively insensitive to GWs at 
low frequencies: in the limit of equal arms, each pair of 
parallel arrows represents the difference of two symmetric 
measurements y.y (t) and y^ (t) that share the same times 
of pulse emission and reception. Taylor-expanding the 
hj?(. ■■',t) terms of Eq. around t and around either 
p\ or p l 2 , we find that yij(t) — yji{t) <x L 3 h'-j. By con- 
trast, the differences of head-to-tail double arrows that 
appear in X sum up to L 2 h'^. Considering monochro- 
matic GWs of frequency Jgw, we see that the GW re- 
sponse is smaller for £ than for X by a factor 7 2tt/q^L 
(~ 0.1 for / GW = 1CT 3 Hz, ~ 0.01 for 10~ 4 Hz). Since 
the response to the LISA secondary noises is approxi- 
mately the same for £ and X [as can be seen using Eqs. 
(j2U)l , discussed below] , ( turns out to be relatively insen- 
sitive to GW. 



5 Most interestingly, a GW-insensitive observable allows the ob- 
servational distinction of a stochastic GW background from in- 
strumental noise |22| . 

6 The interpretation of U as a synthesized observable was already 
clear to F. B. Estabrook (unpublished note). 

7 For unequal interferometer arms, f ex L(AL)(i"^, so the ratio 
between the £ and X responses is ~ AL/L ~ 0.01. 
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Sagnac (a) \ Unequal-Arm Michelson (X) Beacon (P) 




six-link observables \ eight-link observables 



FIG. 4: Mapping first-generation TDI observables into closed arrow loops. The Sagnac (a), unequal-arm Michelson (X), and 
Relay (U) observables of first-generation TDI have a straightforward interpretation as synthesized two-beam interferometers. 
More interestingly, the Beacon (P) and Monitor (E) observables can be seen as four-beam interferometers, with the four beams 
combining into different pairs at the events of initial emission and final reception. To interpret the symmetrized Sagnac (£) 
observables as six-beam interferometers, three different beam pairings must be invoked to explain the cancellation of laser noise 
at emission and reception, and the relative insensitivity of £ to GWs (see main text). 



D. The observables of second-generation TDI 



This interpretation of TDI observables as 2iV-beam 
synthesized interferometers is intriguing, but also trou- 
bling, since it casts a suspicion of arbitrariness on the 
selection of a standard set of observables, and it com- 
plicates exploring the space of all possible combinations. 
Fortunately, the application of the tools of modern alge- 
bra to TDI showed that all first-generation observables 
can be obtained as algebraic combinations of four gen- 
erators 0- This approach was extended to modified 
TDI observables [SlM ITo |. Hi ) , which cancel laser noise 
in rotating LISA geometries, where the Sagnac effect [2l| 
introduces a distinction between light-travel times in the 
two directions. (The Michelson-, U-, P-, and E-type ob- 
servables of first-generation TDI are bona fide modified 
TDI observables, if written with the correct primed and 
unprimcd delay indices; by contrast, the Sagnac observ- 
ables of modified TDI are different, and twice as long as 
those of first-generation TDI.) 

However, the algebraic approach cannot be extended 
easily to the observables of second generation TDI, 
which cancel laser noise in LISA geometries with time- 
dependent armlengths. 8 As pointed out by Cornish and 
Hellings 0, in this situation it is necessary to keep track 
of the order of retardations: for instance, the unequal- 
arm Michelson combination of Fig. 0] would translate to 



With second-generation TDI, the cancellation occurs up to (and 
including) terms proportional to L;L m , which is more than suf- 
ficient for realistic LISA spacecraft orbits. 



yi3'2;322'(£) + 2/231;22'(*) + 2/123;2'(i) + 2/32'l(l) 
- V23l(t) ~ J/l3'2;3 W - 2/32'l;3'3 (0 ~ 2/l23;2'3'3(l), 

where, using the semicolon notation of Ref. |Toj . 

yur,rS) =yuj(t- L ri ), 

Vily^rx (t) — Uilj (t ~ L ri — L r2 (t — L ri )j , 

and so on: the nominal time t is delayed incrementally 
starting from the rightmost delay index n . (A similar no- 
tation with commas instead of semicolons is used when 
the armlengths are constant and the order of the retar- 
dations is not important.) 

Inserting the laser noises Ci in Eq. ©, we see that 
they cancel in pairs, except for the terms from the tails 
of the two paths, 

Cl ; 3'322' — Ci ; 22'3'3(£) — C\ (t) [t ; 3'322' — £;22'3'3jj (H) 

Taylor-expanding the retardations to first order and 
keeping only linear terms in the Li(t), we get 

Ci (i) [(L 2 , + L 2 )(L 3 + L 3 ,)- (L 3 + L 3 , ) (i 2 , +L 2 )], (12) 

where all the Lj and Li are implicitly evaluated at time 
t. (More generally, each retardation index r.; generates a 
residual term proportional to L ri L r for each index rj to 
its left.) 

In short, much like what happened with the simple 
Michelson combination [Eq. @] for unequal-arm geome- 
tries, a laser-noise residual appears in Eq. (|12fl because 
the light-travel times built up along the two interfer- 
ing paths are different; graphically, the tails of the two 
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paths do not match precisely. This is because, although 
both paths contain the same set of links, they do so in 
different orders, and the retardations do not commute 
when the armlengths are time dependent. (This is also 
the reason why the algebraic approach becomes ardu- 
ous for second-generation TDI, where it involves the so- 
lution of polynomial equations for noncommuting vari- 
ables.) As in the upgrade from equal-arm to unequal-arm 
(first-generation) Michelson observables, one solution is 
to compose the two paths so that each goes through each 
arm twice, in different orders [Io|. The residual of the re- 
sulting 16- link combination vanishes up to the first Tay- 
lor order and to the first degree in Li (henceforth, to 
first order/degree). This second-generation unequal-arm 
Michelson observables (known as X\) may be written in 
our notation 

2/l3'2;322'22'3'3(£) + ?/231;22'22'3'3 (t) 
+ 2/l23;2'22'3'3(i) + Z/32'l;22'3'3 (*) 

+ 2/123;2'3'3(*) + 2/32'l;3'3(*) + 2/l3'2;3(*) + 2/231 (*) ^ 
~ Z/32'l(*) - J/123;2' (*) ~ 2/231;22' (*) ~ 2/13'2;322' (*) 

— J/231;3'322'(*) — 2/l3'2;33'322' (*) 

— Z/32'l;3'33'322'(*) ~ 2/l23;2'3'33'322' {t), 

which is related to the X\ defined in Refs. 0, ^| by a 
change of sign and by the use of the opposite convention 
for primed and unprimed indices. 

Second-generation generalizations of all first- 
generation TDI observables were described by Tinto 
and colleagues 0, [n|. For the analogs of the a, (3, 
7, and £ observables (which are formally identical to 
their modified TDI counterparts 0, 0], except for the 
interpretation of the delay indices as noncommuting), 
laser-noise cancellation is not complete, even to first 
order/degree: however, the residuals consist of symmet- 
ric sums of LiL m terms that turn out to be small for 
realistic LISA orbits. 



E. The combinatorial enumeration of TDI 
observables 

Our geometric approach to TDI makes it possible to 
enumerate all the second-generation TDI observables of 
given length. The key to this is the 

Feynman Wheeler 9 Geometric TDI principle: any 



9 "I received a telephone call one day at the graduate college at 
Princeton from Professor Wheeler, in which he said, 'Feynman, I 
know why all electrons have the same charge and the same mass.' 
'Why?' 'Because, they are all the same electron!' [. . .] I did not 
take the idea that all the electrons were the same one from him 
as seriously as I took the observation that positrons could simply 
be represented as electrons going from the future to the past in a 
back section of their world lines." [R. P. Feynman, "The Devel- 
opment of the Space-Time View of Quantum Electrodynamics," 
Nobel Lecture, Dec 11, 1965.] 



2Af-beam Geometric-TDI closed loop can be seen as a 
single beam that travels forward and backward in time 
to meet itself back at its origin. 

For instance, the two-beam equal-arm-Michelson com- 
bination of Fig. 13 (left panel), can be interpreted as a 
single beam that, departs at the initial time t — L3 — Ly. 
travels forward in time to be measured at time t — £3, 
and again travels forward in time to be measured (and 
interfere against itself!) at time t; the beam then moves 
backward in time to be emitted at time t — L21, and 
again moves backward in time to be emitted at the orig- 
inal time t — L,2> — £2 (equal to t — L3 — Ly , since the 
armlengths are equal). This closes the loop, and cancels 
laser noise at all junctions (when we translate graphs to 
formulas, we must remember to give minus signs to all 
the backward-time arrows, drawn dashed in Fig. QJ. 

Once we have established that all n-link loops can be 
represented as a single loop, we can enumerate them 
combinatorially by choosing a starting spacecraft and, 
for n times over, choosing the future or past time di- 
rection, and the leftward (clockwise) or rightward (coun- 
terclockwise) movement direction, in all possible com- 
binations. Each loop can be denoted by the index of 
the initial spacecraft, followed by a string of "L" or "R" 
crested by "— >" for forward-time arrows and by "<— " for 
backward-time arrows; this notation is translated easily 
into strings of link indices crested by their time directions 
(henceforth, link strings). For instance, we would write 

(l)LRLR = ¥3^2 and (l)LRRLRLLR = 3'322'33'2'2 for 
the loops in the left and right panels of Fig. |3 respec- 
tively. 

Not all 3 x 2 2n strings with n links correspond to laser- 
noise-canceling combinations, because the total light- 
travel time accumulated across the loop must be zero 
(for second-generation TDI, zero to first order/degree). 
However, it is quite straightforward to set simple closure 
criteria that identify the true TDI combinations: 

Pre-TDI interferometry. For equal-arm geometries, 

the loop must end at the initial spacecraft (#[L, R] — 

#[R, L] mod 3 = 0, where # denotes the number of oc- 
currences of a symbol in the string), and it must have a 

null total light-travel time (#[L, R] = #[L, R]). We 
denote the combinations that satisfy this property as 
closed. 

First-generation TDI. For unequal-arm geometries 
with generic Li = Ly , the loop must end at the ini- 
tial spacecraft, and satisfy #[/,/'] = (for 
/ = 1,2,3), which yields a null total light-travel time. 
We denote the combinations that satisfy this property as 
\L\- closed. 

Modified TDI. For unequal-arm geometries with 
generic Li =/= Ly, the loop must end at the initial space- 
craft and satisfy #[T] = #[T] (for I = 1, 1', 2, 2', 3, 3'), 
which yields a null total light-travel time. We denote the 
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combinations that satisfy this property as L- closed. 

Second-generation TDI. For unequal-arm geometries 
with generic, time-dependent Li(t) ^= Lp(t), first- 
order/degree laser- noise cancellation is obtained for loops 

that are L- closed, and in addition satisfy #[ I rh, I rh] = 

#[ I rh, I rh] (with I, m — 1, 1', 2, 2', 3, 3'), where a pair 

Im is counted for each link / with itself and with all 

the links m and m to its right, and for each link I with 
all the links rft and m to its right. This condition yields 
a total light-travel time that is null to first order/degree 

(see App.[OJ. For instance, for 3'32'2 we count ¥3' , 3^3, 

3^ 2 , 33, 3 2', 3 2 , and ^ hence the counting 
does not satisfy the property given above. We denote the 
combinations that satisfy this property as L- closed. 

The closure criteria induce useful symmetry properties 
for the link strings: 

Null bigrams. The bigrams L L , R R, L L , and R R 

(or equivalently I I and I I ) always represent combi- 
nations of two yij measurements that sum up to exactly 
zero: any n-link string that contains such bigrams repre- 
sents a combination of length smaller than n. 

Cyclic string shift. Shifting a string cyclically pro- 
duces a combination that has the same closure properties 
as the original, and that differs only by an overall time 
advancement or retardation, and by an advancement or 
retardation applied selectively to the shifted terms, but 
considered negligible at that closure level (for instance, 
for an L-closed loop the additional selective retardation 
would be of leading order/degree L/L m ; for an L-closed 
loop, of a higher order/degree). 

Time and direction reversal. Simultaneously swap- 
ping the primedness and time direction of all link indices 
produces a combination that has the same closure prop- 
erties as the original, and that differs from the original 
by its handedness. 

Cyclic index shift. Shifting link indices cyclically 
(1 — > 2, 2 — > 3, 3 — > 1) produces a combination that has 
the same closure properties as the original, and that dif- 
fers only by a relabeling of indices. Noncyclic index per- 
mutations, on the other hand, produce illegal strings (i.e., 
unconnected loops) unless all indices change in primed- 
ness, in time direction, or in a combination of the two. 10 

String reversal. Reversing a string while swapping all 
"— ►" and "<— " symbols (i.e., reversing all time directions) 
produces a combination that has the same closure prop- 
erties as the original, and that differs only in sign and by 



10 The reflection-like symmetry considered below Eq. 1361 consists 
in exchanging two indices and swapping their primedness, while 
reversing the time direction of the third index. 



an overall time advancement or retardation considered 
negligible at that closure level. 

Splicing. Inserting any link string at any compatible 
point within a link string yields another legal link string. 
Here compatibility means that the spacecraft visited by 
the first loop at the insertion point (for instance, after 1 , 

spacecraft 2; after 1 , spacecraft 3) must be the same as 
the initial spacecraft of the inserted loop. The resulting 
string has at least the closure properties shared by the 
spliced fragments. 

For instance, the unequal-arm Michelson (L-closed) 
string 



i3'23i232'i323'i2'32i 



(14) 



(where the subscripts show the spacecraft visited by the 
loop before and after traversing each link) can be spliced 
at its center with its own L-closed reversal, 



i3'23i232'[i232'i3'23i2'32i323'i]323'i2'32i, (15) 

yielding an L-closed loop that is in fact the second- 
generation TDI Michelson observable. 

At this point it is also useful to give a rule to translate 
link strings into combinations: 

1. Starting at the left end of the string, write 
a yij measurement for each index, according 
to the replacement rules {1, 2, 3, 1', 2', 3'} = 
{V32, 2/13, 2/21,2/23, 2/31,2/12}; attribute a plus sign for 
"— ►" links and a minus sign for " links; 

2. while doing this, build the delay sequence to be ap- 
plied to each new 2/y, adding (from the left) an ad- 
vancement index f before translating each r , and 
a retardation index s after having translated each 



For instance, the 3'322'33'2'2 string would translate to 

^13'2;3 7 (*) + fasijiEFw 
+ ^123;233' W + 2/32'l;2'233' (*) 
— f231;2'233'(*) ~~ 2/l3'2;32'233' W 
~ y32'l;3'32'233' W ~ yi23;2'3'32'233' W- 



(16) 



To accommodate advancement indices, the semicolon no- 
tation of Eq. I|10|) is extended by the advancement rule 



yuj;a(t) = yaj(t + r a ), 



(17) 



with r a defined as the time experienced by light prop- 
agating along link a for emission at time t, and given 
implicitly in terms of L a (t) by r a (t) = L a (t + r a (t)). Re- 
tardations and advancements are applied incrementally 
starting from the rightmost index. 

Equation (|16|l appears more complicated than Eq. ijjjj , 
but it encodes essentially the same y^ combination: the 
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two equations are related by a time retardation, as can 
be seen by evaluating Eq. (fTCj) at the time ^31322' ■ This 
adds 3'322' to the right of the delay sequence for each 
Uij, and since adjacent pairs 11 and 11 cancel in delay se- 
quences (by the very definition of Li and r{), Eq. i|16|) 
then turns into Eq. 10. Conversely, Eq. (|13[l may be 
obtained by applying the translation rule to the string 

3'322'22'3 / 32 / 233'33'2'2, and evaluating the resulting ex- 
pression at the time i ; 33'2'22'233' ■ A slightly more com- 
plicated version (see App. IA 2|) of the rule given above 
yields combinations that are closer to standard TDI 
notation. 



noises is given by 

Vif = Vij ~ 2pmj - pin-., + pm i; i forunprimed l(i, j), 

Vij = Vij ~ P m *j ~ P m j for primed I (i, j) , 

(20) 

where y° p is the optical-path noise in the ytj phase 
measurement, and pnij and pm* are the velocity noises 
of the two proof masses aboard spacecraft i. Because 
prrii-i — prrii and pm* t — pm* have the same time sig- 
nature as laser phase noises, they are canceled in TDI 
observables; thus, all retardations can be removed from 
the unprimed-Z(i, j) expression of y^\ casting it to the 
same form as its primed-i(i, j) counterpart. 



F. Extension of Geometric TDI to six-laser LISA 
configurations 

The extension of our discussion to LISA configurations 
with six lasers introduces three additional laser phase 
noise variables C*(t), C|(i), and C^(t), corresponding in 
Figs.|2and|nito the optical benches on the right-hand side 
of the spacecraft (if we look toward the center) . Equation 
© changes accordingly: 

Vij = Ctl ~ Cj + GWs forunprimed 
Vij = Ci-j — C* + GWs for primed 

here unprimcd and primed link indices l(i,j) correspond 
to c.c.w. and c.w. j/y measurement directions, respec- 
tively. As shown in Ref. [lt| . all arguments and deriva- 
tions valid with three lasers can be applied to a six-laser 
configuration by replacing all with y^' (ry,j in Ref. 
11]) defined by 

Vij = Vij - \( z i]d - z i]d) forunprimed 

Vij = Vij + \( z \i\j - z v) forprimed j), 

where the z%j (t) are the back-plane measurements com- 
paring the phases of the two lasers on each spacecraft, 
and where the notation |Z| removes a prime from I, if 
present. Because % (i) = ±[C*(i) — Cj(t)] [for unprimed 
and primed respectively], the C* disappear from 

Eq. (IBfr. and Eq. @ is restored for the y^' . This justi- 
fies all the developments reported in this paper also for 
six-laser LISA configurations. 

It should be mentioned in this context that the phase 
noise from the random motion of the optical benches en- 
ters the with the same time signature as the laser 
phase noises, and is therefore also canceled by TDI. The 
LISA sensitivity to GWs is then set by the remaining 
secondary noises. Adopting the schematization of the 
measurement process used in most of the TDI literature, 
and the notation used to describe the Synthetic LISA 
simulator [20|. the response of the to the secondary 



III. A GEOMETRIC SURVEY OF 
SECOND-GENERATION TDI OBSERVABLES 

I have written a computer program to list all the 
second-generation TDI observables consisting of 24 or 
fewer j/y measurements. For each even length n, this 
was achieved by enumerating all 2 2 ™ possible LR strings, 
and checking each of them for L-closure, according to 
the counting rule given in Sec. Ill El Already for 24-link 
strings, the combinatorial space is huge, and an exhaus- 
tive search required more than 10,000 CPU hours. The 
resulting list of observables was then reduced to a min- 
imal set by removing all the quasi-duplicates that differ 
only by a sign or by a cyclic string shift. I have kept 
as distinct the observables that differ by a cyclic index 
shift (in first-generation TDI, this would correspond to 
counting X, Y, and Z as separate observables). The 
reduced list of TDI observables is available at the web- 
page www. vallis . org/tdi, annotated with their tempo- 
ral footprint (see Sec. lIIIB|l . number of beams, type, and 
splicing composition (see Sees. IIII Al and lIIIC|l . 

My results are tallied in Table [Q Here the X, U, E, 
and P types represent generalizations of the 8-link ob- 
servables of the same name: X-type (Michelson) observ- 
ables use two arms in both directions, [/-type observables 
use four oriented arms in a relay configuration, i?-type 
and P-type use four oriented arms in beacon and mon- 
itor configurations; the observables tallied under other 
use either five or six oriented arms. Here are the high- 
lights of the survey, which are discussed in more detail in 
the following sections. 

• I find that the shortest second-generation TDI ob- 
servable has length 16. By contrast, modified-TDI 
observables begin at length 8. 

• I recover all 11 the known 16-link second-generation 
TDI observables, previously obtained by Tinto and 



11 In fact, our [/-type combinations do not include the U± observ- 
able given in Ref. Il ll , which achieves laser-noise cancellation at 
the approximate time t — 4L through the sum of four distinct 
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links 


comb. 


space 


unique obs. 


X-type 


P, -E-typc 


U-type 


other 


2-bcam 


4-beam 


6-beam 


8-beam 


10-bcam and higher 


16 


4 


< 10 


48 


12 


18 


18 





3 


27 





18 





18 


7 


< 10 10 


192 





12 


12 


168 


6 


24 


102 


60 





20 


1 


< 10 12 


660 


24 


18 


18 


600 


12 


114 


276 


90 


168 


22 


2 


< 10 13 


2412 





36 


36 


2340 


30 


264 


732 


792 


594 


24 


3 


< 10 14 


12585 


144 


90 


90 


12261 


99 


945 


2676 


4566 


4299 



TABLE I: Geometric survey of second-generation TDI observables. Quasi-duplicates that differ only by a sign or by a time 
advancement are not counted. The X, U, E, and P types represent generalizations of the 8-link observables of the same 
name: X-type (Michelson) observables use two arms in both directions, [/-type observables use four oriented arms in a relay 
configuration, i?-type and P-type use four oriented arms in beacon and monitor configurations; the observables reported under 
other use either five or six oriented arms. The number of beams corresponds to the number of contiguous substrings with the 
same time direction, minimized with respect to cyclic string shifts. All the observables tallied in this table are available in 
extenso at the webpage www.vallis.org/tdi 



colleagues by applying commutator-like de- 

lay operators to the 8-link observables of modi- 
fied TDI. From a geometric viewpoint, all the 16- 
link observables can be understood as self-splicings 
of the 8-link observables of the same type. This 
shows that the former reduce to finite time differ- 
ences of the latter, up to time shifts of the first 
order/degree. It follows that the second-generation 
TDI observables have the same sensitivity of the 
modified TDI observables of the same type, 12 not 
only in the equal-arm limit, but unconditionally. 

See Sec. lIII Al for details. 

• In addition, I obtain alternative forms of the known 
16-link observables. The alternative forms use a 
larger number of beams (e.g., four beams for X, as 
opposed to the standard two), or a different allo- 
cation of links in the beams (e.g., 5 + 4 + 3 + 4 or 

6 + 4 + 2 + 4 for U, as opposed to the standard 

7 + 4 + 1 + 4 structure). The alternative forms, 
too, can be understood as self-splicings of the 8- 
link modified TDI observables of the same kind. 

The alternative forms have the same sensitivity to 
GW signals as the original forms in idealized mea- 
surement conditions, but they can improve on them 
when realistic aspects (such as quantization of the 
phasemeter output and technical noises) are taken 
into account. In addition, the alternative forms 
have a reduced temporal footprint (the difference 
between the times of the earliest and latest phase 
measurements involved in their construction); this 
feature can be advantageous in the presence of gaps 
or glitches in the data, because it reduces the 
extent of defect propagation to the TDI time series. 



yij measurements; by contrast, in Geometric TDI laser noise is 
always canceled by construction bet ween pairs of phase mea- 
surements. However, the Ui of Ref. 1 1 ll has almost the same 
temporal structure as our 1'12'3' 1~2'3'1':L3'2'1'11'3'2' (keep in 
mind that the primedness of our indices is the opposite of Ref. 

ED)- 

Neglecting of course the fact that the modified-TDI observables 
would not cancel laser noise in a flexing LISA. 



See Sees. IIII Al and IIII Bl for details. Appendix IA 41 
gives explicit algebraic expressions for all the 16- 
link observables in terms of the measurements. 

• Second-generation TDI observables are found in in- 
creasing numbers at lengths 18, 20, 22, and 24. A 
minority are of the X , U, E, or P types, while most 
use either five or six oriented arms. 

All 18-to-24-long observables can be understood 
as splicings of modified-TDI observables of length 
8 to 18, sometimes with the inclusion of null Di- 
grams; most, but not all, are self-splicings. I con- 
jecture that all second-generation TDI observables 
of any length can be generated as splicings of two 
modified-TDI observables. 

See Sec. imTil for details. 

• Up to length 24, I do not find any L-closed observ- 
ables of the £ type (defined as having suppressed, 
but nonzero, GW response at low frequencies). I 
conjecture that the £ type is incompatible with 
L-closure. This does not exclude the existence of 
non-L-closed C _ type observables (such as the £i , £2 , 
and C3 described by Tinto and colleagues ^l|) that 
do not cancel laser noise to first order/degree, but 
bring it sufficiently below the LISA secondary noise 
to be useful in practice. 



A. Structure and sensitivity of the 16-link 
observables 

The standard second-generation TDI X observable is 



A+ 6 ' 2 : 3'322'22'3'32'233'33'2'2 (2 beams), (21) 

which is related to the Xi defined in Eq. (|13fl by AT 16 ' 2 = 
A^i- 33'2'22'233" The alternative forms found in our geo- 
metric survey are 



x i6,4,+i . 3'322'33'2'233'22 / 3'32 / 2 

* , _>< m (4 beams), (22) 

x i6,4,-i . 3 '32'233'22'3 / 322 / 33'2'2 
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which differ between themselves only by handedness, and 



X™'*' : 22'22'33'2'23'33'32'233' (4 beams), (23) 

which turns out to have vanishing response in the equal- 
arm limit to both noise and GWs, at all frequencies. To 
see that X\ and the xj 6 ' 4 ^ 1 have all the same GW sensi- 
tivity as the 8-link modificd-TDI X (neglecting of course 
the fact that X would not cancel laser noise in a flexing 
LISA) , we reason as follows. 

As we have learned in Eq. (|15|) , X\ can be interpreted 
as a self-splicing of X with its reversal. If we take X to 
be defined by Eq. 0, we see that 



3 322 33 2 2 — ~^^'2'233' ' 



(24) 



22'3'32'233' = -X. w 



;33'2'2- 



Since the time at the splicing point in 3'322'33'2'2 is 



t. 



2'2:i:S 



7, from Eq. (fT5|l we see that 



3'322'22'3'32'233 / 33 / 2 / 2 ~ X. w 



2'233' 



X. 



;33'2'22'233' ' 



(25) 



here the symbol "~" denotes equality up to selective de- 
lays or advancements of order LiL m in the y^, not spec- 
ified by the formal delay strings. [In this case, these 



spurious delays appear because 22'3'32'233' is only L- 

closed, so the time at the beginning and at the end of 

the inserted string is different by terms of order 

>< 

consequently, the last four y.y observables of 3'322'33'2'2 
are really evaluated at the time £-3'322' 33'2'22'233' > n °t J us * 



2'233' 



Rewriting Eq. I|25(l in terms of X\ and reabsorbing the 
time advancements by evaluating the equation at time 

£;3'322'22'3'3, we find 



Xi ~ x. 



22'3'3 



X. 



(26) 



Thus, up to delays of first order/degree, self-splicings 
produce finite differences of observables. (Indeed, it is 
a well-known fact in the literature on second-generation 
TDI observables that the standard 16-link observables 
are approximately equal to finite differences of the stan- 
dard 8-link observables of the same type.) 

Now, because the individual yy measurements respond 
linearly 13 to GWs and to all instrumental noise sources, 
the strain sensitivity of X to monochromatic sources of 
frequency / at a given sky position is proportional to 



X n (f)/X GWs (f), where X n (f) is the (square-root) spec- 
tral density of noise in X, and X GWs is the Fourier trans- 
form of the GW response function. The constant of pro- 
portionality is SNR/v / 2dtiI, where SNR is the fiducial 
signal-to-noise ratio at which the sensitivity is defined, 
and y/Tobg is the time duration of the observation. Com- 
bining the Fourier-transform time-shifting property with 
Eq. (|26[) . and considering that first order/degree terms 
can be neglected for secondary noises and GWs (which 
are much weaker than the laser phase noises), we see that 
X\ must have the same sensitivity as X: 



( e 2«/At _ l)X n (f) 



X n (f) 



X GWs (f) (e 2 "/ A * - l)X GWs (/) X GWs (/) 



(27) 



with At = t ; 22'3'3 — t. This is true generically for any 
spacecraft geometry, and not just in the equal-arm limit. 

To see that xj 6 ' 4 ' +1 and Xj 16 ' 4 '" 1 , too, have the same 
sensitivity to GWs as X (in the limit of perfect laser 
phase noise cancellation) , it is then sufficient to show that 
they are self-splicings of X and of its (cyclically shifted) 
reversal: 



x WA,+i . 3'322'[33'2'233 / |22']3'32 / 2, 



X 



16,4,-1 



3'3[2'233'|22'3'3]22'33'2'2. 



(28) 



Moving on to the U type, we see (for instance) that 
the three second-generation £/-type observables that use 
the oriented arms 1, 1', 2, and 3, 



3211'231[23|1'1321'1]1', 



3211'23[1'1321'123]11', 



(29) 



3211'[132|1'1231']2311', 



are generated by the self-splicings of the modified TDI 
[/-type observable that uses the same oriented arms, 

3211'231T. Thus, the modified TDI and second- 
generation TDI [/-type observables have the same GW 
sensitivity (in the limit of perfect laser phase noise can- 
cellation). Similar arguments hold for the P- and i?-type 
observables. 

Altogether, we find empirically that all the 16-link 
second-generation TDI observables can be generated as 
self-splicings of the 8-link modified TDI observables of 
the same kind. Conversely, it can be proved that all self- 
splicings of L-closed observables are L-closed (see App. 

rot . 



Linearity is always assumed in the model of measurement used to 
derive TDI. Significant nonlinearity in the phase measurements 
has been explored little, but would probably prove very detri- 
mental to the delicate cancellation of laser phase noise achieved 
by TDI. 



B. Advantages of the alternative 16-link 
observables 

As mentioned above, the alternative forms of the 16- 
link observables can have a smaller time footprint than 
the standard forms. For instance, the standard X\{t) of 
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FIG. 5: TDI response of the second-generation X-type vari- 
ables to the fundamental secondary noises, according to Eq. 
pit , assuming equal arms and proof-mass optical-path noise 
spectral densities given by S? m = 2.5 x 10 _48 [//Hz]_ 
and S° p = LSxlO" 37 [f/Hz] 2 Hz" 1 (following Ref. 
top panel shows the noise response of the standard 16-link X\ 
observable, as compared to the response of the alternative 16- 
link forms of Eq. (I22L which have half as many nulls. The 
middle panel shows the noise response of the two-beam, 20- 
link X\ observables of Eq. 13011 . The bottom panel shows 
the three new noise responses found for 24-link Xi observ- 
ables (solid and dashed: self- and nonself-splicings of 12-link 
observables; dash-dotted: self-splicings of 8-link observables, 
with inclusions). To avoid visual clutter, two of the 24-link 
X\ curves are not plotted beyond / = 1/(21/) (shaded region). 



~ 8L, and each single y^ measurement appears in X\ at 
times displaced by as much as ~ 6L. Thus, X\ will be 
unavailable during the first and the last ~ 136 s (i.e., 
8L) within each LISA data-taking period. Moreover, a 
data gap in a single measurements will appear in the 
X\ time series at four distinct times spanning ~ 6L. By 
contrast, the alternative forms x\ & ' 4,±1 involve 16 y%j 
measurements taken within a time interval of span ~ 6L, 
and the single appear at times displaced by at most 
~ 4L. The gain is significant, if not dramatic. 

The alternative forms for the U-, E- and P-type vari- 
ables also yield footprints gains with respect to their stan- 
dard forms (from ~ 7L to ~ 5L for U, and from ~ 5L 
to ~ 4L for E and P). These gains are possible because 
the alternative observables are obtained, loosely speak- 
ing, by folding the standard versions in time, using both 
time advancements and retardations, as opposed to re- 
tardations only, to arrange the y^ measurements so that 
laser phase noise is canceled at all emission and reception 
events. 

Another advantage of the alternative forms is that they 
can yield an improvement in GW sensitivity in realis- 
tic measurement conditions. In the top panel of Fig. 
13 the dashed curve shows the power spectral density 
(PSD) of secondary noise for the standard X\ observ- 
able, drawn in the limit of equal armlengths. Following 
Ref. [2(j, we assume that secondary noise consists en- 
tirely of proof-mass noise (idealized as stationary and 
Gaussian, with PSD Sf 11 = 2.5 x 10" 48 [//Hz]- 2 Hz" 1 ) 
and of optical-path noise (also stationary and Gaussian, 
with PSD S° p = 1.8 x 1(T 37 [f/Hz] 2 Hz" 1 ). As dis- 
cussed around Eq. 1(771) . the sensitivity to GWs (aver- 
aged over noise realizations) is computed by dividing 
the secondary-noise rms power by the GW transfer func- 
tion. 15 See, for instance, Ref. [6j for plots of the GW 
sensitivities of the first-generation X-, U-, P-, and E- 
type observables, common also (as discussed above) to 
the second-generation observables of the same type. 

Although the noise PSD has nulls at multiples of the 
inverse armlength light-travel time [for the standard X± , 
the nulls are at / = k/(4L)], the sensitivity to GWs 
remains finite in idealized conditions, because the GW 
transfer function displays zeros of the same order at the 
same frequencies. In reality, we should expect a degra- 
dation of sensitivity at these frequencies, because noise 
as a whole can only drop to the level of uncanceled laser 
noise, 16 or of other technical noises (such as quantiza- 



Eq. (|13f) involves 16 j/y measurements taken 14 within the 
interval [min(t 3 '322'22'3'3, i ; 22'3'33'322')i t], for a time span 



14 Although the nominal times of the yu are displaced only by 
~ 7L, we must remember that in a six-laser LISA configuration 
some of the yi^ contain one additional delay, as given by Eq. 

(Hp. 



Roughly speaking, the GW transfer function for the standard 
X\ is found by inserting Eq. J5J into Eq. 1131 and Fourier- 
transforming, assuming a monochromatic source at a fixed sky 
location. The resulting transfer function is usually integrated 
over sky locations. 

However, the simulations of Ref. L'.'fl seem to suggest that un- 
canceled laser noise, under the assumption of perfectly linear 
responses, will display the same nulls as the fundamental sec- 
ondary noises. 
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tion noise), effectively filling in the nulls. The nulls are 
unwelcome also because they imply that a very large dy- 
namic range is needed for sensitive measurement at those 
frequencies. The alternative forms x\ & ' i,±l improve on 
this situation, because their noise PSD and GW transfer 
functions 17 have half as many nulls [at / = fc/(2L)], as 
shown by the solid curve in the top panel of Fig. [SJ 

Similar gains are found for the other 16-link second- 
generation TDI observables. The standard [/-type ob- 
servable of Ref. 1 which we may represent as 18 

321lT32i'123l r 2311', has nulls at / = fc/(3L), while 

the alternative form 3211'231231'1321'H' has nulls only 
at / = k/L. The distribution of the nulls improves also 
for all the alternative forms of the E- and P-type observ- 
ables. 



C. Longer observables 

Although the size of the combinatorial space of LR 
strings scales as 2 2 ™ with increasing string length n, we 
find empirically that the number of valid L-closed combi- 
nations grows only as 2 n (at least up to n = 24); by con- 
trast, the number of the less constrained L-closed com- 
binations grows roughly as 2 3 "/ 2 (at least up to n = 16). 

The majority of the longer second-generation observ- 
ables use cither five or six LISA oriented arms, and there- 
fore do not belong to any of the A, U, E, or P types. 
We do however find new forms for these. Throughout 
this section, we shall consider examples of the X type: 
findings and conclusions are similar for the other types. 

At length 20, we find two A-type observables with two 
beams, 



3'322'22'[22'3'32'233'12'233'2'2, 

- (30) 

3'322'3'3[22'3'32'233']33 / 33'2'2, 

which can be interpreted as the self-splicing of the 8-link 

modified TD I observable 3'322'33'2'2 with its reversal 
22'3'32'233', after the insertion of the double null bigrams 
Tl'Tl and $333' [shown in bold in Eq. Because 
the yij measurements corresponding to the null bigrams 
sum up to zero to first order/degree, we can repeat the 
arguments of Sec. IIII Al to see that the 20-link X\ ob- 
servables of Eq. (|3U|I have the same sensitivity as the 
standard 16-link form. Moreover, reasoning along the 
lines of App. IA 31 we can prove that the self-splicings of 
L-closed strings are L-closed even with the inclusion of 
null bigrams. 



17 In the equal-arm limit, S£ rig ' 16 = 4cos 2 (27r/L)S£ lt ' 16 . 

18 As discussed in note 1111 the observables are not identical: how- 
ever, they do have the same secondary noise PSD and GW trans- 
fer function. 



At length 20 we find also alternative Xi observables 
with four beams, 

3'322 / 22'[33~'|22'3'32'2]2'233'2'2, 

3'322'3'3[3'32'233' 1 2?1 33'33'2'2 , 

L< _L< (31) 

3'33'3[22'3'32'233 / ]33'22'33'2'2, 

3'322'33' 22'[22'3'32'233']2'22'2 , 

and with six beams, 

3'322'33' 22^ [3^3 1 22'3'32'2l2'22'2, 

>< — ' n_ J _><_ (32) 

3'322'33'33'[33'|22'3'32'2]3'32'2, 



all of which are again self-splicings of 3'322'33'2'2 with 
two null-bigram inclusions. The secondary-noise PSD of 
the two-beam, 20-link Ai observables of Eq. I|30|) is shown 
in the middle panel of Fig. and has nulls at / = fc/(6L). 
This PSD is not found for any of the 16-link X\ observ- 
ables. The four-beam observables have the secondary- 
noise PSD of the standard 16-link X\ (the dashed curve 
in the top panel of Fig. EJ, while the six-beam observ- 
ables have the PSD of the alternative 16-link X\ (the 
solid curve in the top panel of Fig. |SJ- 

We find more A-type observables at length 24. Some 
of these, such as 

3'322'22'22'[22'3'32'233']2'22'233'2'2, (33) 

are self-splicings of the 8-link modified TDI X with 
quadruple null-bigram inclusions, and have the same GW 
sensitivity as the standard 16-link X\, They can have 
the same secondary-noise PSDs as the 16- and 20-link 
Ai observables, or a new PSD (the dash-dotted curve in 
the bottom panel of Fig. with nulls at / = fc/(4L) and 
at frequencies given by third-degree algebraic numbers. 
Other 24-link A-type observables, such as 

3'322'22'[22'22'3'32'22'233'133'2'22'2, 

►< (34) 

3'33'322'[22'3'33'32'233'33']33'33'2'2, 

are self-splicings of the 12-link modified TDI A-type ob- 
servables 3'322'22'33'2'22'2 and 3'33'322'33'33'2'2. The 
24-link observables of Eq. I)34|) can be shown to have the 
same GW sensitivity as the standard 16-link Ai in the 
equal- armlength limit. To see this, we parse the 8- and 
12-link A-type variables as self-splicings of the pre-TDI 

(simply closed) observable 3'32'2, respectively without 
and with double-null-bigram inclusions: 

3'3[22'33']2'2, 

3'322'[22'33']2'22'2, (35) 

3'33'3[22'33']33'2'2. 

If the armlengths are equal, the yij measurements from 
the double bigrams shown in bold sum up to zero, and 



give no contribution to the GW and secondary-noise re- 
sponses; if the armlengths are different, the presence of 
the bigrams introduces a preferred direction that dis- 
tinguishes the sensitivities of the 8-link and 12-link ob- 
servables. In the equal armlength limit, the observables 
of Eq. (|34|1 can have the same secondary-noise PSD as 
the standard 16-link X\ observables, or they can have 
two new PSDs (shown as solid and dotted curves in the 
bottom panel of Fig. |3J) with nulls at / = k/(2L) and 
f = (2k+ 1)/(4L), and at / = fc/(4L) and / = fc/(6L), 
respectively. 

At length 24 we find also some X-type observables de- 
rived from noTi-se/f-splicings of 12-link modified TDI X- 
type observables, such as 

3'322'22'[22'3'33 / 32'233'33 / l33'2'22'2, 

m (36) 

3'33'322'[22'22'3'32'22'233']33'33'2'2. 

These non-self-splicings occur between 12-link L-closed 
link strings that differ by a reflection-like, noncyclic index 
shift (e.g., 2 <-> 3', 3 «-> 2', with 1 and 1' reversing time 
direction, if they are present in the string). In the equal- 
armlength limit, these observables have again the same 
sensitivity as the standard 16-link X\, and they have the 
same secondary-noise PSDs as either the standard 16-link 
X\ observables or the self-splicings of 12-link observables 
discussed above. 

We now direct our consideration back to the full set 
of L-closed link strings of lengths 16 to 24. By exhaus- 
tive exploration, we were able to show that all of them 
can be obtained as splicings of two shorter L-closed link 
strings. An explicit specification of the splicings is given 
in the lists of TDI observables available at the webpage 
|www . vail is . org/tdi| A good portion of the L-closed 
link strings are obtained from self-splicings. The non- 
self-splicings occur between L-closed strings that are of- 
ten related by evident symmetries (such as index shifts) . 
Interestingly, these symmetries may occur between in- 
closed strings of different length, as it happens for the 
20-link observable 

3'3[3'|122~'T3'312'213]22 : '33'2'2, (37) 

which can be interpreted as the splicing of the 8-link X 
>< 

string 3'322'33'2'2 with the reversal of the 12-link string 

3'3122'133'T2 7 21. The latter can be obtained from the 
former by reversing the time direction of all 3'3 and 22' 
bigrams, and inserting 1 and 1 bridges as needed to 
keep the arrows consecutive, as shown in Fig. |SJ 

Also interesting is that some splicings result in L-closed 
strings shorter than the sum of the lengths of their com- 
ponents, as in the case of the 22-link observable 

322'21[l'l32i'l23l~3 2'23l'2'3122'l' = 

— >«-L m (38) 

322'211'1321'122'231'2'3122'1', 

which can be interpreted as the splicing of the 16- 
link string 322'2132'231'2>3122T with the 8-link [/-type 
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FIG. 6: The time-reversal plus bridge symmetry is exploited 
by certain L-closed non-self-splicings at length 20 and above, 
such as the link string of Eq. 11371 . 

m 

string 1'1321'123, occurring in such a way that the final 
3 of the latter is erased by the subsequent 3 from the 
former [both bold in Eq. 

We conjecture that all L-closed strings can be obtained 
as splicings of L-closed strings. A possible proof could fo- 
cus on the monadic L-closed strings (those that cannot be 
obtained as the splicing of two L-closed strings), which 
appear at every even length: consider for instance the 

string 3'322'22' • • • ■ ■ ■ 2'2. To prove our con- 

jecture, one would show that monadic L-closed strings 
are never L-closed, as we know to be the case up to length 
24. 



IV. CONCLUSIONS 

I have described Geometric TDI, a powerful new ap- 
proach to understand TDI intuitively, and to interpret 
all of its observables as the measurements of virtual syn- 
thesized interferometers, extending the original intuition 
of Tinto and Armstrong p| and the graphical explana- 
tions of Shaddock ^(| an d Summers |14|- In Geometric 
TDI, observables consisting of n one-way phase measure- 
ments are represented by link strings of length n, and 
can be enumerated exhaustively by listing all possible 
link strings, and then applying simple rules to deter- 
mine which strings correspond to first-generation (|L|- 
closed), modified (L-closed), or second-generation (in- 
closed) observables. A study of the closure symmetries of 
link strings provides clues to the general rule that modi- 
fied TDI observables combine (i.e., splice) into second- 
generation observables, maintaining in most cases the 
same GW sensitivity. 

In addition to its pedagogical value, Geometric TDI 
has the practical interest of providing a systematic 
method to explore the space of second-generation TDI 
observables; such a method was unavailable prior to this 
work. Possible applications include the optimization of 
GW sensitivity (in analogy to the work of Prince and 
colleagues in Ref. 01), and the development of targeted 
noise diagnostics. 

In Sec. IIIII I have used Geometry TDI to survey 
all second-generation TDI observables of lengths up 
to 24. My results (available in full at the webpage 
www.vallis.org/tdi) show that the garden of TDI ob- 
servables contains a wealth of previously unknown spec- 
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imens; this richness only increases with string length. 
In Sec. 1111 Bl I have pointed out how certain alternative 
forms of the standard observables [such as the four-beam 
X\ observables of Eq. I|22l) ] have improved GW sensi- 
tivity in realistic measurement conditions (because they 
have fewer noise and GW response nulls), and reduced 
susceptibility to gaps and glitches (because they have a 
smaller temporal footprint). 

The new observables become possible in Geometric 
TDI because the time advancements and retardations 
of one-way phase measurements are put on the same 
footing, while only retardations were considered in tra- 
ditional TDI. This generalization does not constitute a 
third-generation TDI, but its combinatorial power and 
attractive symmetry justify its addition to the canon 
of TDI. There seem to be no particular challenges to 
implementing time-symmetric observables, especially in 
the framework of post-processed TDI [24[ , whether per- 
formed onboard or on the ground. 

Future work on Geometric TDI should endeavor to 
prove the two open conjectures formulated in this paper: 
namely, that all second-generation TDI observables of 
any length can be generated as splicings of two modified- 
TDI observables, and that there are no second-generation 
TDI observables of the perfect £ type. (It would also be 
interesting to explore a notion of relaxed L-closure that 
includes the Ci , C2 , and £3 variables described by Tinto 
and colleagues 11].) Other promising directions of re- 
search include the optimization of GW sensitivity (can it 
be characterized geometrically?) and the exploration of 
generative rules other than splicing for second-generation 
TDI observables (what is the geometric counterpart and 
generalization of the algebra of observables studied in 
Refs. 0?). 
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APPENDIX A: RULES AND PROOFS 

1. Counting rule for L-closed link strings 

There is nothing magic about this rule, which is nec- 
essary and sufficient to have a null total light-travel time 
at first order/degree for generic armlengths L/(t). To see 



why, set t = at the beginning of the link string; step- 
ping through the string, increment t by an advancement 
Fi (t) (defined as the time experienced by light propagat- 
ing along link I for emission at time t) for each symbol 

I , and decrement t by a delay Li(t) for each symbol 

I . Each advancement F; or delay Li enters the time ar- 
guments of all subsequent -T m 's and L^'s in the string, 
generating terms L{\ x {r m ,—L m } in the first- 

order/degree Taylor expansion of the total light-travel 
time. Since r t {t) = Lj(t+ !"{(*)) ~ L t {t) + Li(t)Li(t), we 
can replace all F;'s with Li's (while keeping the overall 
sum the same to first order/degree) by adding a term 

L1L1 for each I in the string. The counting rule given 
in the main text is then equivalent to requiring that all 
the LiL m terms cancel one by one to yield a null total 
light-travel time. 

2. Translation rule from link strings to 
quasi-standard TDI expressions 

1. Start at a bigram of type I m (there must be at 
least one in every closed loop), shifting the string 
cyclically to move the bigram close to the middle; 

— > 

2. move to the left, starting with I , and write 
down a yij measurement for each index, accord- 
ing to the replacement rules {1, 2, 3, 1', 2', 3'} = 
{V32, 2/13, 2/21,2/23, 2/31,2/12}; use a plus sign for "-»" 
links and a minus sign for "<— " links; 

3. while doing this, build the delay sequence to be ap- 
plied to each new y^, adding a retardation ;r after 
having translated each r , and an advancement ;s 
before translating each *F; 

4. after reaching the left end of the string, go back to 
the link index m in the initial bigram, and move to 
the right, using the same replacement rules; 

5. while doing this, build the delay sequence from 
scratch, adding a retardation ;s after having trans- 
lated each V, and an advancement ;r before trans- 
lating each r . 

Shifting the string cyclically by a few positions can reduce 
the length of the delay sequence. 

3. Proof that all self-splicings of i-closed 
observables are L-closed 

A few lemmas are needed for this proof. Looking back 
to the counting of Irh pairs outlined in Sec. Ill El we de- 
note as prod [string] the polynomial obtained by regard- 
ing the 36 possible Irh pairs as monomials, and summing 

them with coefficients #[Zm, lm] — #[ I rh, I m]. A string 
is L-closed iff prod [string] = 0. 
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Lemma 1. — The prod of an L-closed string is un- 
changed after a cyclic string shift. Consider the shift of 
a single index m from one end of the string to the other. 
Because the string is L-closed, and must therefore have 
#[ I ] — #[ I ], the contribution to prod [string] from the 
index m sums down to zero if m carries a — >, or to mm 
if it carries a <— (and therefore does not multiply itself). 
After the shift, the contribution to prod [shifted string] 
from the index m is the same. Hence this lemma. 

Lemma 2. — The prods of an L-closed string and its re- 
versal are opposite numbers. This is established by notic- 
ing that prod[string] + prod[reversal] = [string] x [string] 
(i.e., the result of counting one Im term, with the appro- 
priate sign, for each I with each m in the string, includ- 
ing itself). However, for a closed string, [string] x [string] 
must be zero; this is because in L-closed strings the con- 
tribution from each index sums down to zero, given that 
each index multiplies every other in equal numbers under 
— > and <— . Hence this lemma. 

Lemma 3. — The cross prod [string] x [reversal] of a 
shifted string and its shifted reversal (i.e., the result of 
counting one Im term, with the appropriate sign, for each 
I in the string and each m in the reversal) is zero. In fact, 
the cross prod is separately zero for each index in the 
string with all the indices in its reversal, again because 
the reversal (as the original string) is L-closed and has 
#[m] = #[m]. Hence this lemma. 

Lemma 4- — All self-splicings can be brought into a nor- 
mal form given by the concatenation ("|") of the shifted 
string and its shifted reversal, for appropriate shifts. This 
does not change the prod of the self-splicing. 

Proof. — Hence, prod [self-splicing] is given by 

prod[shifted string | shifted reversal] = 

prod[string] + prod[reversal] (Al) 
+ [shifted string] x [shifted reversal] = 0, 

because the first two terms are opposite numbers and the 
third term vanishes. Hence the proof. 



4. Algebraic expressions for all second-generation 
TDI observables of length 16 

In this section I give explicit algebraic expressions for 
all the second-generation TDI observables of length 16, 
as found in my exhaustive survey, modulus the symme- 
tries discussed in Sec. Ill El There is considerable arbi- 
trariness in writing these expressions, corresponding to 
the selection of representative link strings in each equiv- 
alence class, to the convention used in translating strings 
to sums of yij measurements, and to the choice of the 
initial time of evaluation for each observable. Here I list 
link strings in a normal form whereby each string begins 
with the largest continuous substring of forward-time in- 
dices; I adopt the translation rules given in Sec. Ill E| [just 
below Eq. I|15|) ]: and I adjust the time of evaluation to 



minimize the length of the longest delay sequence in the 
expression. 



X type. The standard 16-link observable X 16 ' 2 is ob- 
tained by applying the translation rules to the link string 

S'322'22'3'32'233'33'2'2, and evaluating the resulting ex- 
pression at the initial time i ; 3'322'22'3'3: 



yi32;322'22'3'3 + £/231;22'22'3'3 + 2/l23;2'22'3'3 
+ 2/321;22'3'3 + ?/l23;2'3'3 + 2/321;3'3 + 2/l32;3 + 2/231 

— J/32'l — 2/l23;2' ~ 2/231:22' ~ 2/l32;322' ~ 2/231;3'322' 

— 2/l32;33'322' _ 2/321;3'33'322' — 2/l23;2'3'33'322' \ 



(A2) 



the alternative 16-link observables x 16 ' 4 ^ 1 can be writ- 
ten from the lin k stri ngs 3'322'3'32'233'2?33'2'2 and 
22'3'322'33'2'23 7 32'233', evaluated at times t 
and i22'3'322'33 7 ' respectively: 



3'322'3'32'2 



2/l32;322'3'32 7 2 + y231;22'3'32 7 2 J'l23;2'3'32 7 2 
+ y321;3'32 7 2 ^132;32 7 2 + V231^3 ~ V^l-^Pl — 2/l23;2 
^ 1/231 ^ 2/132;3 + 2/l23;23'3 + 2/321;2 7 23'3 ~~ V 2?,1^2 ! 2 , 3' 
~ ^132;32 7 23'3 — 2/321:3'32 7 23'3 — 2/l23;2'3'32 7 23'3' 



(A3) 



yi23;2'3'322'33 7 ~"~ 2/321 ; 3'322'33 7 2/i32 ; 322'33 7 
+ ^231;22'33 7 + 2/l23;2'33 T ' %21;33 7 — V2?,l^i¥ ~ 2/l32;3 7 
~ 2/32'l ~ 2/l23;2' + J/l32;3 7 22' + 2/231;33 7 22' ~ 2/321;33 7 22' 
_ yi23;2'33 7 22' — y231;22'33 7 22' — yi32;322'33 7 22' ' 



(A4) 

last, the null 16-link observable xj 6 ' 4,0 can be written 
from the link string 22'22 ; '33'2'23'33'32'233', evaluated 

at time ^22'22'33 7 2 7 2 : 

^123;2'22'33'2'2 + y321;22'33'2'2 + ^123;2'33'2'2 
y321:33'2'2 — f231;33'2'2 — 2/l32;3'2'2 — 2/321;2 7 2 — 2/l23;2 
yi32;3 7 y23l{W 2/i32;3'33' + 2/231;33'33' — y321;33'33' 
— yi23;2'33'33' ~~ 2/231;22'33'33' — 2/l32;322'33'33' ' 

(A5) 

Expressions for the A-type observables that use the arms 
1, 2 and 1, 3 are obtained by cyclic index shifts. 

U type. The three 16-link [/-type observables that use 
the oriented arms 1, 1', 2, and 3 correspond to the link 
strings 



3211'231231'1321'H', 
3211'23l'132i'12311' , 



(A6) 



3211'1321'1231'2311'; 
applying the translation rules and evaluating at the times 
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^32ii'23i2> ^3211'231'1j an< ^ ^32ii'i32iJ respectively, yields 



2/231;211'2312 2/l23;ll'2312 + 2/312;l'2312 
+ 2/213;2312 — 2/l23;2312 — %31;3l2 — 2/3l2;T2 — 2/l23;2 

— 2/231 + 2/213;F3 ?/312;TF3 + 2/231;311'3 ~"~ 2/l23;2311'3 

— 2/213;2311'3 — f 312:1'2311'3 — 2/213;ll'2311'3 ' 



(A7) 



2/231;211'231'l + 2/l23;ll'231'l + 2/312;l'231'l 
+ J/213;23X'l ~ fl23;231'l ~ 2/231;31'l + 2/213;! + 2/312 
+ 2/231;3 ' 2/l23;23 — 2/213;23 — 2/312;l'23 — 2/l23;ll'23 
— 2/231;211'23 — 2/312;3211'23 — 2/213; 13211 ' 2"3~ ' 

2/231;211'132F + 2/i23 ; ll'l32F %12;1'132F 
+ 2/213;132~ + 2/3l2;32F 2/231;2F + 2/i23;F — 2/213;F 
- 2/312 - 2/l23;l - 2/231;21 + 2/ 2 13;F321 ~ 2/l23;F321 
_ 2/231;2F321 — 2/312;32F321 — %13;132F321 ' 



(A8) 



(A9) 



the third expression is closest to the Ui given in Ref. ^l| 
(but see note II H) . Expressions for the [/-type observ- 
ables that other sets of oriented arms (i.e., {2, 2', 1,3}, 
{3, 3', 1,2}, {1,1', 2', 3'}, {2,2',1', 3'}, and {3, 3', 1', 2'}) 
are obtained by cyclic and noncylic index shifts. 

E type. The three 16-link E-type observables that use 
the oriented arms 1, 1', 2', and 3 correspond to the link 
strings 



H'132'1'12' 3 1'2'311'3 2', 
1'11'2'311'T?132 7 1 7 1? < 3 , 
n 7 ?3"l'l32'l'l2^3iT'"3 2 7 ; 



(A10) 



applying the translation rules and evaluating at the 

times ^n/132'1'12' > ^1'11'2'3TF3> an< ^ ^H'2'31'132 7 ' res P ec " 

tively, yields: 

2/312;l'132'l'12' 2/213;132'l'12' 2/312;32'l'12' 

2/231:2' 1' 12' — 2/321;2'l'12' — 2/213;Fl2' ~ 2/312;T2' + 2/32'l 

— 2/231 + 2/213;F3 2/321;2'l'3 _ 2/231;2'l'3 ~~ 2/312;32'l'3 

— 2/213;132'l'3 2/231;31'132'l'3 — 2/321;31'132'l'3' 

(All) 

2/213;ll'2'311'3 + %12;1'2'311'3 2/213;2'311'3 
+ 2/321;311'3 — 2/231;311'3 — 2/3l2;TF3 — 2/213;F3 + J^ 231 

— 2/32'l + 2/312;T2' + 2/231;312' — 2/321;3T2' — 2/213;2'3T2' 

— 2/312;l'2'3T2' 2/321;2 7 ll'2'3T2' — 2/231;2 7 ll'2'3T2' ' 

(A12) 



2/312;l'2'31'132 7 + 2/213;2'31' 132 7 + 2/321;31'132 7 
~ y231;31'13F + J^lS^lF 2/3i2;3F 2/231;F — 2/321;F 
— 2/21'3 ~ 2/312:1' + 2/321;2 7 ll' — 2/231;2 7 ll' — 2/312;32 7 ll' 
_ 2/213; 132 7 !!' + 2/231;31'132 7 ll' — 2/321;31' 132 7 !!' ' 

(A13) 

the third expression is closest to the Ei given in Ref. Hi!. 
Expressions for the other possible sets of oriented arms 
(i.e., {2, 2', 1, 3'} and {3, 3', 1', 2}) are obtained by cyclic 
index shifts. 

P type. The three 16-link P-type observables that use 
the oriented arms 1, 1', 2, and 3' correspond to the link 
strings 



3'1'H' 2 3'11'3'213' 2 1'12, 

2ii'i^^i'i23'i : 'T?n 7 3', 



(A14) 



applying the translation rules and evaluating at the 

times ^3/1/11/23'TF' ' : 2ii'i3 7 2FT' an d ^211'23'i'is 7 ' respec- 
tively, yields 



2/l32;l'll'23'Fi 7 1" 2/213;ll'23'Fi 7 ' 2/312;l'23'Fl 7 
2/213;23'TF — 2/l23;23'TF + 2/l32;TF — 2/3i2 ; TF ~ 2/213;F 

— 2/l3'2 + 2/l23;23' + 2/312;l23' — 2/l32;l23' + 2/l23;23'T23' 

— 2/213;23'123' — 2/312;l'23'T23' — 2/l23;ll'23'l23' ' 

(A15) 

2/l23;ll'13 7 2FT + 2/3l2;l'13 7 2FT + 2/213;13 7 2FT 
+ 2/312;3 7 2FT — 2/l32;3 7 2FT 2/l23;FT — y213;FT — 2/312;T 

— 2/123 + 2/l32;3 7 2 + 2/213;l'3'2 — 2/l23;l'3'2 2/l32;3 7 21'3'2 

— 2/312;3 7 21'3'2 — 2/213:13 7 21'3'2 — 2/l32:l'13 7 21'3'2 ' 

(A16) 

yi23;ll'23'l'13 7 2/312; 1'23' 1' 13 7 ' 2/213;23' 1' 13 7 

— 2/l23;23'l'13 7 2/l32;l'13 7 2/213; 13 7 ' 2/3l2;3 7 — 2/l32;3 7 
+ 2/l23;2 — 2/213;2 — 2/312;l'2 — 2/l23;ll'2 + 2/l32;3 7 211'2 

_ 2/312;3 7 211'2 _ 2/213;13 7 211'2 — 2/l32;l'13 7 211'2' 

(A17) 

the third expression is closest to the Pi given in Ref. hll. 
Expressions for the other possible sets of oriented arms 
(i.e., {2, 2', 1', 3} and {3, 3', 1,2'}) are obtained by cyclic 
index shifts. 
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